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Abstract 

The thermalization of a longitudinally expanding color glass condensate with Bjorken boost 
invariant geometry is investigated within microscopical parton cascade BAMPS. Our main focus 
lies on the detailed comparison of thermalization, observed in BAMPS with that suggested in 
the "Bottom-Up" scenario. We demonstrate that the tremendous production of soft gluons via 
gg — > ggg, which is shown in the "Bottom-Up" picture as the dominant process during the early 
preequilibration, will not occur in heavy ion collisions at RHIC and LHC energies, because the 
back reaction ggg — > gg hinders the absolute particle multiplication. Moreover, contrary to the 
"Bottom-Up" scenario, soft and hard gluons thermalize at the same time. The time scale of thermal 
equilibration in BAMPS calculations is of order a~ 2 {\.na s )~ 2 Q~ 1 . After this time the gluon system 
exhibits nearly hydrodynamical behavior. The shear viscosity to entropy density ratio has a weak 
dependence on Q s and lies close to the lower bound of the AdS/CFT conjecture. 
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I. INTRODUCTION 



It is of great interest to understand the mechanisms of quick equlibration of quark matter 

^. , 

the large values of the elliptic flow v 2 measured in Au+Au collisions [2J, [3[ . When the quark 
gluon matter becomes sufficiently dilute due to the strong longitudinal expansion pQCD, 
bremsstrahlung processes are essential for momentum isotropization 4], 0]. 

The importance of pQCD bremsstrahlung was first raised in the "Bottom-Up" scenario 
which describes the thermal equilibration of a Color Glass Condensate (CGC) [3, Q] 
characterized by a saturation scale Q s . The main idea of the "Bottom-Up" scenario is that 
while the hard gluons with transverse momenta of order Q s degrade as the condensate evolves 
in space time, soft gluons with transverse momenta much smaller than Q s are populated 
due to pQCD gg — > ggg bremsstrahlung and start to dominate the total gluon number. As 
soon as the radiated soft gluons achieve thermal equilibrium and build up a thermal bath, 
the hard ones begin to loose their energy to the thermal bath and subsequently thermalize. 
A parametric time scale for overall thermalization is given by r th ~ a~ 

Thermal equilibration of gluonic matter with an idealistic CGC initial condition in pres- 
ence of pQCD gg <-> ggg processes is investigated in the present work for the first time 



within transport calculation using the parton cascade BAMPS 



12j. In this work we apply 



a one dimensional expanding geometry, which is adequate to the assumption of Bjorken 
boost invariance [3]. For the initial gluon distribution of CGC we employ an idealized and 
boost-invariant form 0, Q f(x,p) = —jjr-8(y — rj)Q(Q 2 s —pr), with N c = 3 for SU(3) and 
c ~ 0.4 {14I flo) ]. The initial gluons are produced at eigentime r ~ This form of initial 
condition is a simplified form of CGC inspired by earlier studies in jol, [h]] . Particles with 
transverse momenta p± > Q s are missing here. Their presence would modify the thermal- 
ization process, but as well lead to an even faster thermalization of the hard momentum 
region. Including a more realistic distribution in the hard region of transverse spectrum is 
a task for future work. 
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FIG. 1: ^(central rj bm),a s = 0.3 



FIG. 2: Scaled effective temperature 



II. RESULTS: THERMALIZATION OF A CGC 



The way that thermalization proceeds within the parton cascade calculations does not 
resemble the "Bottom- Up" scenario jf]]. The strong parametric enhancement of the total 
gluon number at early times, as predicted by "Bottom-Up" scenario, is not observed in the 
cascade calculations. Instead, gluon annihilation occurs during the first 0.3 — 0.75 fm/c for 
Q s = 2 — 4 GeV, as shown in Fig. [TJ which is due to domination of the 3 — > 2 processes 
at the early stage of evolution. This indicates that the initial CGC is oversaturated for 
the chosen values of a s and Q s . The gluon number begins to increase after t ~ 0.3 — 0.75 
fm/c (depending on the value of Q s ) when the system is close to kinetic equilibrium and a 
quasi-hydrodynamical cooling sets in. 

After the expansion starts, energy flows immediately into both the soft (p t < p so ft = 1-5 
GeV) and hard momentum region (p t > Q s ) where the populations rapidly increase, as 
demonstrated in Figj3](a) for calculations with a s = 0.3 and Q s = 3 GeV. The presence of 
a thermal bath of soft gluons seems not to be a necessary condition for the equilibration of 
hard gluons, since they achieve an exponential shape on a short time scale and almost as 
quick as the soft gluons, as Figj3|a) demonstrates. Again, this is different from the picture 
invoked in the "Bottom-Up" scenario. 

At t = 0.5 fm/c the transverse spectrum looks almost thermal. After t w 0.5 fm/c the 
cooling of the system sets in and the spectra evolve in a way, characteristic for a hydrody- 
namical expansion, as shown in Figj3](b) for a s = 0.3, Q s = 3 GeV. 

We extract the time scale when the system is more or less thermalized from Fig. [2} where 
the scaled effective temperature Tt 1 ^ 3 is shown. We define the time scale of thermalization 
as the time, at which Tt 1 ^ 3 becomes a constant, which corresponds to a one-dimensional 
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FIG. 3: Transverse momentum spectra at early times (a) FIG. 4: Ratio of the shear viscosity 
and for later times (b) to the entropy density. 

ideal hydrodynamical expansion. The times extracted from Fig. [2]( constant a s = 0.3) read: 
1.2, 0.75, 0.55 fm/c (resp. for Q s = 2,3,4 GeV). For fixed Q s = 3 GeV and various a s 
the thermalization times are extracted in same manner: 1.75, 1.0, 0.75 fm/c (resp. for 
a s = 0.1,0.2,0.3). 

From the values presented here we observe that the time scale of thermalization is of 
order a^ 2 (ln a s )~ 2 Q~ 1 . The 1/Q S scaling is consistent with the "Bottom-Up" result. The 
dependence on a s proves to be weaker than the "Bottom-Up" prediction but consistent 
with the scaling obtained in [if]]. The reason for the quick thermalization observed in our 
calculations is the ft 



uon bremsstrahlung, which favors large-angle radiation due to the LPM 
suppression |4|, |5J, |12 ] . Thus, the pQCD gg <-> ggg processes are dominant for the thermal 
equilibration. 

The ideal hydrodynamic behavior of the system at late times indicates that the ratio 
of the shear viscosity to the entropy density is small. Figure H] shows the rj/s ratio in 
the calculations with a s = 0.3 and Q s = 2, 3, and 4, respectively. The shear viscosity is 
calculated using the Navier-Stokes approximation: rj = | (T xx + T yy — 2T ZZ ). The entropy 
density is calculated by s = An — n ln(A) , where n is the gluon density and A = n/n t ^ 
denotes the gluon fugacity. 

The rj/s ratio is nearly constant and has a weak dependence on Q s , r]/s ~ 0.15, which is 
exactly the same as that obtained in full 3+1 dimensional BAMPS calculations with a s = 0.3 

n 

and minijets type initial conditions for Au+Au collisions at RHIC energies [17|. This verifies 
that the rj/s ratio determines the behavior of the late dynamics and, thus, depends only on 
the coupling a s , but not on initial conditions. The smallness of the rj/s ratio, which is close 
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to the AdS/CFT lower bound 18], corresponds to the efficiency of the pQCD gg <-> ggg 
processes in thermal equilibration, because the 7]/s ratio is inversely proportional to the 
total transport collision rate which is for gg ^ ggg processes 6 — 7 times larger than that of 
gg — > gg collisions [if] ]. 

III. CONCLUSION 

Using the parton cascade BAMPS, we found that several aspects of the real thermalization 
might be different compared to the "Bottom-Up" scenario. A strong increase in soft gluon 
number, as predicted in the "Bottom- Up" scenario, is not observed and thermal equilibration 
of soft and hard gluons occurs roughly on the same time scale due to 2 — > 3 and 3^2 
processes, respectively. In agreement with the "Bottom-Up" scenario, the thermalization 
time proves to be proportional to Qj 1 , however, its proportionality to a~ 13 ^ 5 is not seen but 
is much weaker: r t h ~ (a s In a s )~ 2 Q~ l . The differences arise because the back reactions of 
bremsstrahlung, 3 — 2 processes, play a significant role. They are completely absent in the 
"Bottom-Up" scenario. 

The quick thermalization and the smallness of the rj/s ratio observed in the present 



calculations wit 
vious studies 



h t 



he CGC initial conditions are consistent with the findings from the pre- 
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171 | using the Glauber-type minijets initial conditions. This 
demonstrates that independent of the chosen initial conditions, the pQCD bremsstrahlung 
processes (and the back reactions) dominate the dynamical equilibration and then keep the 
system behaving like a nearly perfect fluid. 
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